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TEMPORARY $T0HA$8 0, % Uetttiteatf specified % parameters located ia ^ 

PARAMETERS 3 00 F 00 aP a is the location of the first 

word of the auxiliary routine* 
k 00 F 00 (A.,,)F Storage positions A., A. = A n + i 

contain dependent variables 
y. (i * 0, lj .o., a-l). Initial 
values of y. must "be placed herig 
originally. 

5 00 F 00 (1 *A Q )F Storage positions B^ B ± * A * 

(B n -A^) + i, contain hy! where 
h is the increment. Storage 
positions C,* C. - Au + Sfl^-A ) + i 
are used for temporary storage* 

Originally B. and C. ®ust contain 

• 2 x x 

... • lay .J and h hi*, respectively. 

Locations 7). s !>.. = A^ + 2(l~~A ) + i 
contain values of h^y" found 
hy the auxiliary subroutine. 

6 00 F 00 of n Is the number of differential 

equations to he integrated. 

7 00 F 00 eP * e is the address of the tolerance, 

€« € Is a number such that if, * 
for all i, the ahjsoltfte value 
©f the difference between two 
successive approximations tQ 
hy" is less than 2 -^ e, the 
iteration process is assumed to 
have converged » 



AIS1LIAHY SUBROUTINE 



DURATION 



DESCRIPTION 



The purpose of the auxiliary subroutine is to take the 
quantities y. and hy* (i * 0^ 1, +**, R-l), whose locations 
are determined by the parameters A. ana" B i and compute the 
quantities h 2 y w * The quantity h is the increment of the 
independent variable. The auxiliary subroutine must be 
located in a sequence beginning with location a where a 
is specified by the parameter S3« In integrating over 
one 8tepj the integration subroutine will call in the 
iterations necessary for convergence. 

The time in milliseconds for one time through the subroutine 
is 

Q.<?k + 2.17n + 1(0.25 + 5°8n + D) 
where I - number of iterations for an interval 

D = time of auxiliary routine 

n = number of equations . 
This routine will integrate n second order differential 
equations of the types 

(1) y^f^ y i' " 

i — Of lj « • 9 f B""l 



> y n~l* 



yk> y-I > ■ ° • > yi_i * 



'0* '1 



n-l' 



When entered with the entry 



50 P 

26 integration routine, 

the subroutine will carry out one integration step 
replacing values of y. 



3-1 



and hy. ' 



i , j-l 



by y 



i*j 



and 



hy! 



3 



THE INBEPENDElfC VARIABLE If the independent variable x occurs in the functions 



f , or if it is required during an integration as an 

i 
index, then it may be obtained by integrating the 

equation y" - 0. The independent variable x is then 

treated as an additional dependent variable, for which 

no additional computing is required by the auxiliary 

subroutine. If the independent variable does not 



REMARKS 



-3- 

appear in any of the f . 's, but is merely wanted for 
indication purposes, it is quicker to use a simple 
counter in the main routine. 

The numerical procedure is given by equations 2, 3» 
and k-* 

< 2 > *U* = r itJml + ( h / 2 ) t tf3 ^ + (h/2> f ljJ(k 

if, for all i, the relation 



y. 



IjtJ*^ 



Ijtj 



f = f . 



lj^jfe Ij J,yfc*"i """ 



is satisfied. 



In the above equations? 

i indicates the variable of integration, 

J indicates the step of integration, 
and k indicates the number of iterations. 

Sufficient conditions to guarantee that the procedure will 
converge and that numbers stay in bounds can be stated, 
but use of these conditions leads to extremely short 
intervals. In practice, if units of measure and a 
value for h are chosen so that numbers are within bounds 
within the auxiliary routine the value of h so chosen 
is usually satisfactory for convergence of the numerical 
procedure* 
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